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ABSTRACT 


We present the results of a search for cold gas at high redshift along QSO lines-of-sight carried out without any a priori assumption 
on the neutral atomic-hydrogen content of the absorbers. To do this, we systematically looked for neutral-carbon (Ci) /1/11560,1656 
transition lines in 41 696 low-resolution QSO spectra (1.5 < Zcm < 4.46) from the SDSS-II - Data Release seven - database. Ci 
absorption lines should indeed probe the shielded gas in the neutral interstellar medium of galaxies more efficiently than traditional 
tracers such as neutral atomic-hydrogen (H i) Damped Lyman-o" (DLA) and/or Mg ii systems. 

We built up a sample of 66 CI absorbers with redshifts in the range 1.5 < z < 3.1 and rest-frame equivalent widths 0.1 < Wr(/11560) < 
1.7 A. The completeness limit of our survey is M/.jim(/il560) - 0.4 A. Ci systems stronger than that are more than one hundred-times 
rarer than DLAs at Zabs = 2.5. The number of C i systems per unit redshift is found to increase significantly below z = 2. We suggest 
that the C i absorbers are closely related to the process of star formation and the production of dust in galaxies and that their cosmic 
evolution is driven by the interplay between dust shielding and the evolution of the ultra-violet background at ~ 10 eV. 

We derive the neutral atomic-hydrogen content of the C i systems observable from the southern hemisphere from VLT/UVES spec¬ 
troscopy and find that a majority of them are sub-DLAs with AfHl) ~ 10^** atoms cm“^. The dust content of these absorbers is yet 
significant as seen from the redder optical colours of the corresponding background QSOs and their reddened spectral energy distri¬ 
butions, with £(B-V) values up to ~ 0.3. The overall A(Hl) distribution of Ci systems is relatively flat however. As a consequence, 
among the Ci systems classifying as DLAs there is a probable excess of strong DLAs with log A(Hl) >21 (atoms cm“^) compared 
to systematic DLA surveys. While the dust content of such systems is significant, their dust-to-gas ratio must still be limited. Indeed, 
strong DLAs having large amounts of shielded gas and dust producing stronger reddening and extinction of the background QSOs if 
they exist should have been missed in the current magnitude-limited QSO sample. 

We study empirical relations between VL-fCi), A(Hl), EfB-V) and the strength of the 2175 A extinction feature, the latter being de¬ 
tected in about 30% of the C i absorbers. We show that the 2175 A feature is weak compared to Galactic lines-of-sight exhibiting the 
same amount of reddening. This is probably the consequence of current or past star formation in the vicinity of the C i systems. We 
also find that the strongest C i systems tend to have the largest amounts of dust and that the metallicity of the gas and its molecular 
fraction is likely to be high in a large number of cases. The C i-absorber sample presented here hence provides ideal targets for detailed 
studies of the dust composition and molecular species at high redshift. 
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1. Introduction 

Understanding the mechanisms of star formation at high red- 
shifts is central to our knowledge of how galaxies formed 
and subsequently evolved chemically. This is specially true at 
z ~ 2 when the cosmic star-formation activity was highest 
dMadau & Dickins^l2014h . Stars form out of cold gas, metals 
and dust in molecular clouds (e.g. lSnow & McCallll2006h in the 
interstellar medium (ISM) of galaxies. In turn, the radiative and 
mechanical feedbacks from stars have a strong impact on the 
physical state of the ISM. Studying the ISM at high redshifts, 
and in particular deriving the physical properties of the diffuse 


* Based on data from the Sloan Digital Sky Survey (SDSS) and ded¬ 
icated follow-up observations carried out at the European Southern 
Observatory (ESO) under programmes 082.A-0544 and 083.A-0454 
(RI. Ledoux) using the Ultra-violet and Visual Echelle Spectrograph 
(UVES) installed at the Nasmyth-B focus of the Very Large Telescope 
(VLT), Unit-2 (Kueyen), on Cerro Paranal, Chile. 


molecular phase in galaxies, is therefore crucial for understand¬ 
ing how stars formed in the early Universe. 

The best way to derive physical properties accurately 
is to detect the t racers of the cold gas in absorption (see 
iMuller et al.l l2014l) . The neutral, shielded and possibly cold 
gas clouds at high redshifts can be searched for in the ra¬ 
dio domain by targeting th e neutral atomic-hy drogen (Hi) 21- 
cm absorption line (e.g., iGupta et"^ l2009l) . However, sys¬ 
tematic - blind - surveys have to await the increased sensi- 
tivity of new faci l ities such as Me erKAT/SKA and ASKAP 
(iBooth et alj|2009l: iDuffv et aDl2012l) . On the other hand, such 
a gas can be efficiently traced in the optical wave-bands by 
detec t ing the redshifted H i damped L yman-g (iProchaska et al.l 
2005 : ^terdaeme et al.ll2009bl 2012h and/or s trong Mgii (e.g.. 


Ouider et al.l 201 ll: iBudzynski & Hewetlll20lll) lines imprinted 


in the spectra of bright enough background sources such as 
QSOs or the rapidly fading y-r ay burst (GRB) afterglows (for 
the latter, see lFynbo et al.ll2()(M and references therein). 
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Damped Lyman-o' systems (hereafter DLAs) observed in 
QSO spectra have column densities of A(H i) > 2 x 10^° atoms 
cm“^ and are known to contain most of the neutral gas in the 
Universe in the redshift range 0 < z < 5 (see lWolfe et al.ll2005l 
for a review). It has been shown however that DLAs typically 
prob e warm {T ^ 3000 K) and diffuse (ttn < 1 cm ~^) neutral gas 
(e.g.. iPetitiean et ^120001: ISrianand et al.ll2()T^ . The metallic- 
ity of DL As is generically low, i.e., on an av erage about 1/30* 
of Solar (iPettinj 120061: iRafelski et al.l 1201 2h and their dust-to¬ 
gas ratio is typically les s than one-tenth of what is observed in 
the Galactic ISM (e.g.. IVladilo et al1l2008h . This probably ex¬ 
plains the low detection rates of molecular hydrogen (H 2 ) in 
DLAs where only about 10% of the QSO lines-of-sight inter¬ 
cept H^-be aring gas down to a limit of AfH?) ~ 10*^ molec ules 
cm“^ (e.g.. iNoterdaeme et akll^OSat iBalashev et aI1l2014t for 
searches f or H? in DLAs originating from the ho st galaxies of 
GRBs, see lLedoux et al.OoO^lKruhler et al.ll2013h . 

Based on the observ ed correlation betwe en metallicity and 
dust depletion in DLAs (iLedoux et al.ll2003h . DLAs with high 
metallicity are expected to con tain more dust and t herefore to 
exhibit larger H 2 fractions (see IPetitiean et al1l2006h . However, 
even in DLAs with the highest metallicities typical dust signa¬ 
tures like reddening of the background QSOs, the 2175 A ex¬ 
tinction feature (hereafter also called ultra-violet [UV] bump), 
or diffuse interstellar bands, are not apparent. Even in the rare 
cases with H 2 detections, the inferred molecular fractions are 
low and typical of what is seen in Galactic diffuse atomic gas 
with /(H 2 ) = 2A(H 2)/[2A(H2)-hMH i) 1 ^ 0 01 and often much 
lower than this (see lNoterdaeme et al.ll200^ . The primary rea¬ 
son for this is that the cold and dusty phases are missed proba¬ 
bly because of their reduced cross-s ections relative to that of th e 
more pervasive warm neutral ISM (IZwaan & Prochas^l2006l) . 
Direct evidence for the relatively small physical sizes 0.15 pc) 
of H 2 -detected clouds in DLAs recently came from the observa- 
tion of partial coverage of the QSO broad- line emitting region 
(IBalashev et al.ll^l iMKlimenko et al.ll20L5h . 

H 2 -detected clouds in DLAs are found to have kinetic tem¬ 
peratures in the range T ~ 70 - 200 K and particle densities 
ne ~ 1 - 100 cm“^ (e.g.. ISrianand et aT1l2005l) . When detected, 
H 2 i s usually coincident with neutral atomic carbon (Ci; see 
also lGe & Bechtoidlll999l) . This is due to the fact that the ion¬ 
ization potential of neutral carbon (11.26 eV) is similar to the 
average energy of the Lyman-Werner photons that dissociate H 2 . 
Therefore, shielding of UV photons is essential for these species 
to remain at detectable levels. Carbon monoxide (CO) has long 
escaped detection even in DLAs with detected H 2 , down t o 
N(CO) ~ 10 '^ molecules cm“^ (see, e.g.. IPetitiean et al]|2002l) . 
This is not surprising since CO, with a dissociation energy of 
11.09 eV, needs to be even more shielded than H 2 and Ci to 
be detected. After CO UV absorption bands were detected for 
the first time at high redshift, in a sub-DLA toward s the QSO 
SDSS J 143912.05H-111740.6 (ISrianand et al.ll2008bh . it became 
clear that the best place to detect CO in absorption at high red¬ 
shift are the systems with strong C i absorption. Following this 
strategy allowed us to detect carbon monoxide subsequently i n 
five additional systems (INoterdaeme et al.ll2009aLl201()aLl201 ll) . 

In Galactic translucent interstellar clouds, CO starts to be 
produced in significant amounts when neutral atomic carbon be¬ 
comes the dominant ca rbon species and a larg e fraction of hy¬ 
drogen turns molecular (ISnow & McCallll2006h . The strength of 
the CI absorption is expected to be such that it could be detected 
even in a low resolution spectrum. We therefore embarked in 
a systematic search for C i absorption in QSO spectra from the 
SDSS-II - Data Release seven (hereafter DR 7) - database. In 


this paper, we present the results of this search and the basic 
properties of the detected C i absorbers. Note that here we will 
equally refer to QSO absorption-line systems detected through 
CI absorption as “C i systems” or “C i absorbers”. 

In Sect. |2] we describe our selection and identification of C i 
absorbers. We discuss the properties of the sample in terms of in¬ 
tervening C i-absorber number per unit redshift, proximate sys¬ 
tems, and CI rest-frame equivalent widths, in Sects. l3Jll3.41 and 
13.11 respectively. We then assess the impact of the C i absorbers 
on their respective background QSOs both from the observed 
QSO optical colours (Sect. ITTI) and the reddening these systems 
induce on the QSO spectral energy distributions (Sect. 14.2b . In 
Sect. |5] we present the H i column-density distribution of the C i 
systems from VLT/UVES spectroscopy. In Sect. 0 we investi¬ 
gate empirical relations between neutral atomic-carbon and neu¬ 
tral atomic-hydrogen contents, QSO reddening and the strength 
of possible 2175 A extinction features (whose measurements are 
described in Sect. 14.3b . We summarise our findings and conclude 
in Sect.|7] 

Throughout this paper, we assume a standard A cold dark- 
matter cosmology with Hq - 70 kms“* Mpc^', Qa = 0.7 and 
Om = 0.3. 


2. CI absorption-line selection and identification 


We systematically searched for Ci absorption lines in high- 
redshift QSO sp ectra fr om the Sloan Digital Sky Survey 


(lYork et a l.||2000t) - DR 7 ( Abazaiian et al.ll2009t) - quasar cat- 
alogue dSchneider et al.ll^ioll . This survey imposed an i-band 
magnitude cut of 19.1 for QSO candidates whose colours indi¬ 
cate a probable redshift smaller than ~ 3. The spectra cover the 
wavelength range 3800-9200 A at a resolving power R ~ 2000. 

We implemented a dedicated IDL procedure to detect and 
identify absorption-line features in SDSS QSO spectra auto¬ 
matically. Since SDSS spectra are log-lambda binned, the pix¬ 
els have constant velocity size (» 69 km s '). This makes it 
straightforward to cross-correlate the spectra with an emission- 
or absorption-line template. We used the method introduced 
by INoterdaeme et al.l (1201 Obh to search for [O iii] /1/14959,5007 
emission and Mgii absorption lines. We first normalized the 
spectra iteratively using Savitzky-Golay filtering. This consists 
of smoothing the spectra by convolving them with a Savitsky- 
Golay kernel that preserves the sharp QSO emission-line peaks 
but ignores narrow features such as metal absorption lines, bad 
CCD pixels and sky emission-line residuals. Deviant pixels and 
their neighbours are then masked out and the resulting data is 
convolved again in the same way, and so on and so forth. This 
procedure has the major advantage as no a priori assumption is 
required about the functional form of the QSO continuum (i.e., 
power law or other) and in addition it is computationally ex¬ 
tremely fast. We then cross-correlated the normalised spectra 
with a synthetic profile of C 1/1/11560,1656 absorption lines. We 
looked for the positive correlation signal together with peak ab¬ 
sorptions detected at more than 2cr and 2.5cr, respectively, and 
differing by less than a factor of three. The simultaneous detec¬ 
tions of the Si 11/11526 and Al ii /11670 absorption lines were re¬ 
quired to support the identifications of the two features as C i and 
hence minimize the probability of chance coincidence. Spurious 
detections (~ 50%), most of them close to the detection limit, 
were identified visually and removed from the sample. In total, 
we find 66 systems, one of which is shown in Fig.[T] 

The search for C i lines was limited to the regions of the spec¬ 
tra redwards of the QSO Lyman-a emission line to avoid the 
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spurious coincidences that are frequent in the Lyman-o' forest. 
The wavelength range above 7200 A was also not considered 
in the search to avoid regions of the spectra heavily affected by 
residuals resulting from the sky emission-line subtraction. We 
requested that the search window encompasses the wavelengths 
of the Si 11/11526 and Alii/11670 absorption lines of the puta¬ 
tive CI systems so that the validity of a system does not rely 
solely on the detection of two transitions (see above). For a given 
line-of-sight, the redshift lower bound (Zmin) of the C i search is 
therefore the largest value between z = 3820/1526 - 1 ^ 1.50 
and (1 -H Zem) X 1215/1526 - 1 ^ 0.8 x Zem - 0.2, where Zem 
is the QSO emission redshift. The redshift upper bound (Zmax) 
of the search is the smallest value between Zem + 0.1 (to not ex¬ 
clude a priori proximate systems with infalling velocities of up to 
+5000 kms^') and z = 7200/1656 - 1 ^ 3.35. In order to avoid 
too many false positives at low signal-to-noise (S/N) ratio, we re¬ 
quested the median S/N ratio per pixel to be larger than four for 
a given spectrum to be actually scanned. This resulted in a sam¬ 
ple of 41 696 QSOs with 1.5 < Zem < 4.46 whose spectra were 
searched for intervening or proximate C i absorbers. Note that 
we did not initially reject Broad Absorption-Line (BAL) quasars 
because our procedure follows the QSO continuum locally and 
can detect narrow absorption lines embedded in broad and not- 
fully saturated troughs. Regions of deep absorption are de facto 
avoided when we study the number of intervening C i absorbers 
per unit redshift in Sect. l3.3l as they have low S/N ratio per pixel. 

Table [1] summarises our C i sample which we refer to in the 
following as the overall sample. QSO names with J 2000 coor¬ 
dinates are given for each absorber. No line-of-sight is found to 
feature more than one systerrQ. The SDSS plate and fibre num¬ 
bers as well as the MJD are also provided in the table as use¬ 
ful cross-references. The QSO emission redshifts derived by the 
SDSS team are indicated together with the absorption redshifts 
and rest-frame equivalent widths of the C i lines. The latter were 
carefully determined by us for each individual system. In the 
next column of the table, we specify the average S/N ratios per 
pixel in the regions the two Ci lines are located. At Zabs > 2.2, 
the Lyman-a line of the systems is also covered by the SDSS 
spectra. We therefore provide in the table a determination of the 
total neutral atomic-hydrogen colu mn density of these system s 
following the method developed bv iNoterdaeme et aP (l2009bl) . 
The reliability of the latter is confirmed in a number of cases with 
follow-up high-resolution VLT/UVES spectroscopy (see the last 
column of Table[T] and Sect.|5ll. This column also has additional 
A(H i) measurements for several low-redshift systems for which 
Lyman-a absorption is not covered in the SDSS spectrum. 

3. Sample properties 

It must be noted that in this survey C i systems are found without 
any presumptions on the presence of neutral atomic hydrogen, 
i.e., the C i systems found in this work do not necessarily have to 
be DLAs. Moreover, DLA absorbers can be observed in SDSS 
spectra only when their redshifts are larger than 2.2 while Ci 
lines can be identified down to Zabs - 1 -50. As stated before, the 
reality of the identified C i systems in our sample was checked 
by visual inspection. Therefore, we believe the C i detections are 
secure. We discuss the completeness of the survey in Sect. 13.21 

* Note however that there is a second Ci system towards 
SDSS J 234023.67-005327.1, which, with z^bs = 1-36, falls below the 
redshift cut-off of our survey. This system happens t o be detected in 21- 
cm absorption henc e is also related to cold gas (see lGupta et alll2009l : 
iKanekar et al.ll20Toh . 


3.1. Line equivalent widths 

Because we will rely on them in the analysis, we here seek to ver¬ 
ify the robustness and accuracy of the equivalent-width measure¬ 
ments of C i-absorption lines performed in SDSS spectra. For 
this purpose, we plot in Fig. the measured rest-frame equiv¬ 
alent widths of the Ci/i/ll560,1656 lines versus each other. It 
appears that except in one case the strengths of the two lines are 
in the expected range. This gives confidence in the derived val¬ 
ues and their associated uncertainties. In the case of the outlier 
seen in the lower part of the plot (i.e., at Zabs = 1.526 towards 
SDSS J 125552.60+223424.4), an unrelated blend to the /il656 
line is a probable reason for the observed deviation. 

All the systems are located within about 2cr of the boundaries 
defined by the optically-thin regime on one hand and the relation 
expected for heavily saturated profiles on the other hand. We 
note that because of their large equivalent widths {W^ ^ 0.4 A) 
most of the absorbers, especially those with equivalent-width ra¬ 
tios consistent with the optically-thin regime, are probably made 
of numerous velocity components. 

Here, we assumed that the C i ground-state is solely respon¬ 
sible for the absorption lines while in reality the absorption from 
the two fine-structure energy levels of the neutral-carbon ground 
state (^Pi and ^P 2 ) could in principle contribute mildly to the 
measured equivalent widths. However, this will affect the equiv¬ 
alent widths of both of the 41560 and /il656 lines in the same 
way so that any departure from the assumed relations due to this 
blending will be small. 

3.2. Completeness 

Before discussing the number of C i absorbers per unit redshift 
{tici, see the following section), we first need to estimate the 
completeness of the sample. Given the resolving power R of 
the SDSS spectra, the C i/il560 line rest-frame equivalent width 
limit is given by: 

1560 , 

W,,iin,(/H560) ^ „ X X S/N-' (1) 

K 

where ii - 2 is the number of standard deviations above which 
the peak absorption must be detected and S/N > 4 is the limit 
on the signal-to-noise ratio per pixel at the corresponding line 
position. Note that the FWHM of the lines is sampled by two 
velocity pixels of constant value. Our survey should therefore be 
complete down to Wi jim(/H560) ^ 0.4 A. 

We checked the exact level of completeness of our survey at 
this equivalent-width limit by implementing the following pro¬ 
cedure. For this purpose, we used the same data set that we used 
for the calculation of nci in Sect. 13.31 i.e., the same quasar sam¬ 
ple, the same [Zmin,Zmax] values and the same mean S/N > 4 
limit. We then randomly selected 1000 QSO spectra and in¬ 
troduced an artificial C i system of rest-frame equivalent width 
VTr(/ll560) at different positions in the spectra where the local 
S/N ratio at both C i lines is larger than four. The distribution of 
the equivalent width ratio of the /il560 and /il656 lines is as¬ 
sumed to be a normal distribution with a dispersion correspond¬ 
ing to what is seen in Fig.|2] Note however that neither the equiv¬ 
alent width ratio nor the exact number of artificial systems used 
in the simulation has any significant impact on the completeness 
we infer. 

We implemented about 40 000 C i systems that we sought 
to recover by using the same automatic procedure described in 
Sect. 12 We varied Wr over the range 0.1-1.0 A and defined the 
completeness as the ratio of the number of recovered systems to 
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the total number of systems introduced in the spectra. The results 
are displayed in Fig. [3 It can be seen that the completeness is 
larger than 80% for VTr(/ll560) > 0.4 A. 

3.3. Number of absorbers per unit redshift 

We calculated the sensitivity function, g(z), of our survey, i.e., 
the number of lines-of-sight probing a given redshift z and hav¬ 
ing S/N > 4 at the expected positions of both Ci lines. This 
function is shown in Fig. 01 It combines together the [Zmin,Zmax] 
pairs, previously defined in Sect. |2] for all the lines-of-sight. We 
further excluded the regions with velocities relative to the QSO 
emission redshifts smaller than 5000 km s ', which could, in 
principle, be influenced by the quasar (see Sect. I3.41 i. Note that 
the uncertainties on Zem are of the order of 500 km s“' (see, e.g., 
iParis et al.l2012h and therefore are small enough not to affect the 
statistics. The total statistical absorption path length probed by 
the QSO sample over z = 1.50 - 3.35 is Az ss 13 000 with an 
average redshift (z) = 1.9. From Fig. 0] it is apparent that the 
sensitivity of the survey is an increasing function towards lower 
redshifts. One therefore expects a larger number of absorbers to 
be found at z < 2. This is what is observed in practice as indi¬ 
cated by the redshift histogram of the detected intervening C i 
absorbers over-plotted on the same figure. On the other hand, 
the small number of Ci systems found at Zabs > 2.2 (i.e., eight 
systems out of a total of 66 systems or, equivalently, 12% of the 
sample) is striking. 

To investigate this further, we calculated the number of inter¬ 
vening Ci absorbers per unit redshift, nci, in two redshift bins 
of roughly equal total absorption path length (with boundary 
redshift z = 1.9). Here, we only considered the systems with 
rest-frame equivalent widths above the completeness limit of the 
survey, i.e., Wr(/11560) > 0.4 A. The results are summarised 
in Table |2] and shown in Fig. |5] In the table, we also separated 
the strongest from the weaker absorbers (around a median rest- 
frame equivalent width of 0.64 A) but there is no obvious differ¬ 
ence between the redshift evolution of these two groups. 

The tiQi ~ 1.4 X 10^^ we measure in the higher redshift bin 
(1.9 < Zabs < 3.35), taking into account the effect of incom¬ 
pleteness estimated in Sect. 13.21 implies that C i systems with 
Wr(/il560) > 0.4 A are more than one hundred-t imes rarer than 
DLAs at Zabs = 2.5 (see iNoterdaeme et al.ll2012h . An evolution 
of tici, with nearly thrice as many systems below z = 1.9 than 
above that, is also observed. Compared to the redshift behaviour 
of a non-evolving population, this is significant at the 4.3cr level. 
Such an evolution is interesting and should be studied further as 
it depends on the balance between dust shielding and the UV ra¬ 
diation field. This may imply a strong evolution of the shielding 
of 10 eV photons by dust between z = 2.5 and z = 1.5. 

3.4. Proximate systems 

There are 14 Ci systems with velocities relative to the QSO 
emission redshifts smaller than 5000 kms“'. These could be 
associated with the QSO host galaxy or nearby environment. 
Six systems even have absorption redshifts larger than the cor¬ 
responding QSO emission redshifts (by up to ~ 4000 km s ') 
which is difficult to explain by large peculiar velocities in inter¬ 
vening systems. Imposing the same data-quality cuts and min¬ 
imum equivalent widths as in the previous section, we find 
that the incidence of C i absorbers at small velocity differences 
from the quasars is consistent with that of intervening systems. 


However, associated errors are large due to small number statis¬ 
tics. 

Because of the clustering of galaxies around the massive 
QSO host galaxies, an excess of proximate C i systems could be 
expected. However, C i is a fragile species which can easily be 
photo-ionised by the intense UV radiation emitted by the QSO 
engine. Interestingly, the lack of a significant excess of p roxi- 
mate systems was also observed bv iProchaska et all (l2008h con¬ 
sidering DLAs. In this case, the abundance of proximate DLAs 
is only a factor of two larger than that of the overall DLA pop- 
ul ation. A similarly l ow over-abundance factor was observed 
bv iFinlev et al.l (1201 3h for strong DLAs with logV(Hi) >21.3 
(atoms cm“^). This is much less than what is expected based on 
clustering arguments alone. 

In this work, we do not observe that the properties of proxi¬ 
mate CI systems are different from those of intervening C i sys¬ 
tems. This is true for redshift, Ci equivalent-width. Hi content, 
reddening and UV bump-strength distributions. Nevertheless, 
because of the possibly different origin of these absorbers and 
the possible requirement of strong dust shielding from the 
nearby QSOs, we discriminate in the following proximate Ci 
systems from the rest of the population and comment, whenever 
possible, on proximate C i systems of interest. 

4. Evidence for dust 

4. 1. QSO opticai coiours 

In order to assess the impact of the C i absorbers on their back¬ 
ground QSOs and check for the existence of dust in these sys¬ 
tems, we first consider the observed colours of these QSOs and 
compare them with the colours of the overall QSO population 
used as a control sample. 

In Fig. we show the distributions of (g - r), (r - i) and 
(r-z) colours for the 41 696 QSOs whose spectra were searched 
for CI absorption. In the upper panels of this figure, it is ap¬ 
parent that the lines-of-sight with detected C i absorption do not 
distribute in the same way as the other lines-of-sight. They are 
displaced altogether towards redder optical colours compared to 
the average loci of the QSO redshift sequences. The effect is 
most easily seen in the lower panels of Fig. |6] which compare 
the colour histograms of the two QSO populations (i.e., the Ci- 
detected QSO sample and the overall QSO sample). The two- 
sided Kolmogorov-Smirnov test probability that the two sam¬ 
ples are drawn from the same parent distribution is as small as 
^ lO^'". The typical colour excess is ~ 0.15 mag, i.e., about five 
times lar ger than the mean (r - z) colour excess of 0.03 mag de¬ 
rived by IVladilo et al.l (l2008h in Zabs ~ 2.8 DLAs from SDSS 
DR 5. A similar result for DLAs was found by iKhare et'aP 
(1201 2h based on (g - i) colours of SDSS-DR7 QSOs. This is 
clear evidence for the presence of dust among C i absorbers. 

4.2. Reddening 

Motivated by the unequivocal signature of dust in the form of a 
colour excess of the background QSOs with detected C i absorp¬ 
tion, we now aim at constraining the properties and the nature of 
dust in these systems. 

For each of the 66 QSOs with foreground C i ab¬ 
sorbers, we derived the QSO redden ing, EfB-V), f o llowing 
the same approach a s used i n, e.g. , ISrianand et al.l (l2008ah 
and INoterdaeme et al.l (l2009ai l2010ah . First, we corrected the 
QSO spectra for Galactic reddening using the extinction maps 
from ISchlegel et akl (Il998l) . We then fitted the spectra with the 
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SDSS QSO composite spectrum from IVanden Berk et al.l (1200 ih 
shifted to that QSO emission redshift and reddened with ei¬ 
ther a Small Magellanic Cloud (SMC), Large Magellanic Cloud 
(LMC), LMC2 supe r-shell or Milky Way (MW) extinction law 
( Gordon et al.ll200^ at the Ci-absorber redshift. Our procedure 
is illustrated in the left panel of Fig. |2l The fit with the smallest 

value indicates the most representative extinction law for a 
given absorber. The latter is specified in Col. ’Best fit’ of Table[3] 
and the corresponding ^(B-V) value is given in the preceding 
column. 

For each QSO line-of-sight exhibiting C i absorption, we de¬ 
fined a control sample (hereafter denoted as “C.S.”) made of 
SDSS-DR7 QSOs from the searched sample having an emission 
redshift within +0.05 and a z-band magnitude within +0.1 mag 
from those of the QSO under consideration. In some instances, 
this resulted in a sample of less than 30 QSOs in which case we 
increased the above maximum magnitude difference by steps of 
0.01 until the number of QSOs in the control sample reached (or 
exceeded) 30. We then applied to each QSO spectra from the 
control sample the exact same fitting procedure as described in 
the previous paragraph. Table [3] lists the number of QSOs and 
the median reddening and standard deviation of the distribution 
of ^(B-V) values in each control sample (see the upper right 
panel of Fig. |2]for an illustration). The values given in Table [3] 
correspond to the most representative extinction law previously 
determined for that particular C i-detected QSO spectrum. 

In the left panel of Fig. [8] we show the histogram of red¬ 
dening for the sample of 66 QSO lines-of-sight with detected 
CI absorbers compared to the cumulative control sample (cal¬ 
culated as the sum of the normalized distributions of individual 
control samples). As in Sect. 14. II an offset between the two sam¬ 
ples is apparent. The mean reddening induced by C i systems is 
0.065 mag. A tail in the histogram of the C i-detected lines-of- 
sight is observed, with ^(B-V) values up to ~ 0.3 mag. 

4.3. The 2175 A extinction feature 

A number of C i-detected QSO spectra are best-matched by an 
extinction law exhibiting the absorption feature at rest-frame 
wavelength 2175 A. In order to measure the strength of this 
UV bump (denote d Ahnmn), w e follo wed a prescription similar 
to the one used by iJiang et ^ (1201 Oh where the observed QSO 
spectrum was fitted with the SDSS QSO composite spectrum 
reddened via a parametrized pseudo-extinction law made of a 
smooth component and a Drude component. However, here we 
fixed the wavelengt h and width of the bump to t he Galactic val¬ 
ues determined by iFitznatrick & Mas's^ (l2007h . Both of these 
quantities indeed show little variation from line-of-sight to line- 
of-sight through the Galaxy and the Magellanic Clouds. This 
then limits the number of free parameters and prevents the fit 
from diverging towards very wide and shallow solutions which 
could be non-physical. Indeed, imperfect matching of the ob¬ 
served QSO continuum by the smooth co mponent is expecte d 
due to intrinsic QSO-shape variations ('e.g.. lPitman et al.ll2000l) . 

The fitting process is illustrated in the left panel of Fig. |7] 
The shaded area represents the measure of the bump strength. 
This is the difference between the above best-fit function and the 
same function but considering only its smooth component (i.e., 
with the Drude component set to zero). Abump values are listed 
for each absorber in Table [3] As previously done for the deter¬ 
mination of reddening (see Sect. 14.21) . we also defined a QSO 
control sample whose measured Abump distribution is shown in 
the lower right panel of Fig. |7](i.e., for the given QSO emission 


redshift). Table [3] gives for each control sample the median and 
standard deviation of this distribution. 

The histogram of bump strengths in the C i-absorber sample 
is displayed in the right panel of Fig. [8] One can see from this fig¬ 
ure that more than a quarter of the C i systems feature absorption 
at 2175 A. This strengthens the result from the previous section 
that significant reddening of the background QSOs by dust is in¬ 
duced by some of the C i absorbers. We will come back to this 
and quantify the effect in Sect. |6] 

In the following, we shall use the median ^(B-V) and 
Abump values of the control samples, i.e., (£’(B-V))c.s. and 
(Abump)c.s.j to define the exact colour excess and bump strength 
towards a given Ci-detected QSO line-of-sight: ^(B-V) = ^(B- 
V)measured “ <£'(B-V)>c.s., and likewise for Abump- These zero- 
point corrections are usually almost negligible (see Table |3]l. In 
addition, the standard deviations of ^(B-V) and Abump values in 
each control sample provide an estimat e of the uncertainty due 
to intrinsic QSO-shape variations (see iPitman et al.ll2000l) and 
hence the significance of the reddening induced by each C i ab¬ 
sorber and the significance of associated 2175 A absorption, re¬ 
spectively. 

5. HI content 

As part of a spectroscopic campaign which we will describe in 
a companion paper, we followed up the C i absorbers from the 
overall sample which are observable from the southern hemi¬ 
sphere using VLT/UVES. We present in Fig. |9] the Hi column- 
density distribution of this C i-absorber sub-sample (referred to 
in the following as the H i sub-sample) and compare it with the 
distribution of A(Hi) from systematic DLA and/or sub-DLA 
surveys. 

We secured Hi column-density measurements for most 
of the systems in the overall sample which have a declina¬ 
tion of (5 < +28 deg, i.e., 14 out of 16 systems at redshifts 
Zabs >1.8 (the two exceptions being the lines-of-sight towards 
SDSS J091721.37+015448.1 and J233633.81-105841.5) and 
four out of eight systems at Zabs ~ 1-75 (see Table [B- While 
HI column densities derived from UVES spectroscopy are usu¬ 
ally more accurate, the last two columns of Table [1] show that 
they confirm those derived directly from SDSS spectra as testi¬ 
fied by the five systems at Zabs >2.2 where this measurement 
could be done from both datasets. Eor this reason, we here com¬ 
plement our UVES measurements with the values we derived 
using SDSS spectra for the three systems at Zabs > 2.2 for which 
high-resolution spectroscopic data are not available because the 
background QSOs are too far North for the VET to observe them. 
The HI sub-sample thus comprises a total of 21 systems. 

In Eig.|9] we compare the observed H i column-density distri¬ 
bution of C i-selected abs orbers with that of H i-sel ected DLAs 
(from SDSS DR 7 as welk iNoterdaeme et al.ll2009H) . In this fig¬ 
ure, we also show the expected number of sub-DLAs using the 
fitted distribution function from iProchaska et'BI (l2014ll . We find 
that a large fraction of the Ci absorbers have neutral atomic- 
hydrogen column densities slightly below the conventional DLA 
limit (V(Hi) = 2 X 10^“ atoms cm“^) and therefore classify as 
strong sub-DLAs. However, the fraction of Ci absorbers among 
sub-DLAs is much less than among DLAs indicating that effi¬ 
cient shielding is much more difficult to obtain below the DLA 
limit. Though rare, the existence of C i absorbers with low neu¬ 
tral atomic-hydrogen column densities supports the presence of 
dust in these systems. The dust-to-gas ratio in these systems has 
to be high enough so that the absorption of UV photons by dust 
allows CI to be present in large amounts. 
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No Cl system is found with logN(Hi) < 19 (atoms cm“^). 
This is a regime where shielding of UV photons becomes ex¬ 
tremely difficult even in the presence of dust. However, it is pos¬ 
sible that such systems are missed in our search. Indeed, as seen 
in Fig. [To] there is a trend for neutral atomic-hydrogen column 
density to increase with C i equivalent width. The gradually de¬ 
creasing completeness fraction of the survey below lTi (/ll560) » 
0.4 A (see Sect.|2]i would therefore preclude low A(H i) systems 
from appearing in our sample. 

Above the DLA limit, where the incompleteness fraction 
of our survey is less of an issue, it appears that Ci-selected 
absorbers do not follow the statistics of Hi-selected DLAs. 
Although the number of Ci systems with logA(Hi) > 20.3 
(atoms cm“^) is small (i.e., only 9 systems), it is apparent that the 
overall A(H i) distribution of C i systems is relatively flat. A two- 
sided Kolmogorov-Smirnov test applied to all absorbers with 
log A(H i) > 20.3 (atoms cm“^) gives a probability of only 17% 
that the two distributions come from the same parent population 
(see inset of Fig. 0. This could be explained by a larger num¬ 
ber of velocity components in higher Hi column-density gas, 
thereby increasing the probability of detecting Ci. Moreover, 
large amounts of shielded gas are probably the consequence of 
the line-of-sight passing through the absorbing galaxy at small 
impact parameter, in which case we can expect the A(H i) distri¬ 
bution to be flatter than that of the overall DLA population. 

The strongest DLA found among the Ci absorbers in the 
Hi sub-sample has A(Hi) = 10^''* atoms cm“^. It is however 
located at z^bs ~ ^em- Even ignoring proximate systems, it is 
yet surprising that three intervening DLAs with A(Hi) > 10^' 
atoms cm“^ are present in such a small absorber sample. From 
DLA statistics alone, the probability of randomly selecting three 
DLAs that strong out of a sample of six DLAs is only 6%. There 
is therefore a probable excess of strong DLAs among Ci ab¬ 
sorbers. While the dust content of these systems is significant 
(see Sects. ItTI and 14.21) . their dust-to-gas ratio must be limited. 
Indeed, dust reddening and extinction of the background QSOs 
will inevitably reduce the incidence of strong and dusty DLAs 
in magnitude-limited QSO samples. This implies that the actual 
proportion of strong DLAs among C i systems in general is likely 
to be even higher than what we here found out. 


6. Empirical relations in the sample 

Based on the results presented in the previous sections, we now 
study the existence of empirical relations between the differ¬ 
ent quantities measured in this work: neutral atomic-carbon and 
neutral atomic-hydrogen contents, the reddening C i-selected ab¬ 
sorbers induce on their background QSOs and the strength of 
possible 2175 A extinction features. Because the C1/11560 tran¬ 
sition line is weaker and hence exhibits less saturation than 
CI /11656, we adopt the equivalent width of the former as a proxy 
for the amount of neutral atomic carbon in the systems. 

In Fig. [To] we plot the C1 41560 rest-frame equivalent width 
versus logA(Hi) for the Ci absorbers from the Hi sub-sample. 
Both quantities appear to be weakly correlated. A Kendall rank- 
correlation test indicates the significance of the correlation to be 
1.8cr only. There is therefore a tendency for strong DLAs to have 
larger C i equivalent widths but at the same time, for a given H i 
column density, the C i content can vary substantially from one 
system to another. Large values of VTr(/ll560) are observed in 
DLAs but also in sub-DLAs. The fraction of shielded and prob¬ 
ably cold gas could actually be large in some of these sub-DLAs. 
From the optically-thin approximation applied to the C141560 


absorption line and assuming the ionization equilibrium relation, 
N{C i)/ A(C ii) ~ 0.01, valid for the cold neutral medium (CNM; 
see, e.g.. lJenkins & Tri^l201 ll) . a lower limit on the gas metal- 
licity can be derived: 


[X/H] ^ 18.35-flog 


I Wr(41560) \ 

\0.01 xA(Hi)j 


( 2 ) 


For lTr(41560) = 0.4 A and logA(Hi) = 20 (atoms cm“^), the 
metallicity should be of the order of Solar. More generally, the 
dashed and dashed-dotted curves in Fig. [TOl were calculated us¬ 
ing the above equation assuming metallicities of one-tenth of 
Solar and Solar respectively. Within measurement uncertainties, 
most of the C i systems lie in between these two curves. If the 
medium probed by the line-of-sight is a mixture of cold and 
warm gas, the metallicity of the systems will be even higher. 
However, if part of the hydrogen is in molecular form, the metal¬ 
licity will be lower. For the whole C i-absorber sample, Eq. [2] 
implies a metallicity distribution ranging between [X/H] = -1.4 
and metallicities in excess of Solar, with a median value of 
[X/H] -0.5. This means that the metallicities of Ci absorbers 

would on average be at least ten ti mes larger than those of typ¬ 
ical DLAs (for the latter, see, e.g.. lRafelski et all (l2012h l. This 
should be confirmed by accurate measurements of metal column 
densities. 

In Eig. [TT| we display the relation between ^(B-V) and 
logA(Hi) among the Ci absorbers from the Hi sub-sample. 
Here again, the data points are highly scattered. Most of the sys¬ 
tems are associated with low albeit consistently non-zero QSO 
reddening. Since most of the A(H i) values are relatively low, the 
measured amounts of reddening, with median ^(B-V) ~ 0.045, 
are actually remarkable. This departs from what is observed in 
the overall DL A population where the reddening is us ually negli¬ 
gible (see, e.g.. lVladilo et al.|[2008l : lKhare et al.ll2012h . The latter 
authors have shown that DLAs at Zabs ~ 2.8 typically induce a 
reddening ^(B-V) ~ 5 x 10“^ mag. Apart from a few outliers, 
most of the Ci systems have reddening properties consistent 
with those of the Galactic ISM. This is represented in Eig. [TTIbv 
the solid line. In the Galaxy, the reddening induced along a line- 
of-sight is indeed directly proportional to the neutral atomic- 
hydrogen column density, with EfB-Vl/Al H) = 1.63 x 10“^^ 
mag atoms ' cm^ (iGudennavar et al.ll2012h . Only two Ci sys¬ 
tems with large A(H i) are more consistent with what is seen in 
typical DLAs and/or along SMC lines-of-sight, where the above 
ratio is smaller than in the Galaxy. Two other C i systems with 
low A(Hi) might also deviate from the Galactic relation be¬ 
ing consistent with a ten-times larger ratio. We caution however 
that the uncertainties on the reddening measurements are fairly 
large. If real, this would imply in the latter systems the exis¬ 
tence of a grain chemistry more evolved than in the Galaxy with 
a larg er fraction of big grains over very small grains (e.g., iPeil 
I 1992 I) . This is opposite to the trend observed in the Magellanic 
Clouds. In such systems, strong 2175 A absorption is expected. 
Interestingly, this is what is observed in practice as these two C i 
systems exhibit two of the strongest three Abump values of the H i 
sub-sample. 

To investigate the characteristics of the C i absorbers further, 
we look in Eig.[T2|in more detail at the properties of dust in these 
systems. In the left panel of this figure, ^(B-V) and lTr(41560) 
are found to be correlated with each other at the 4.4cr signif¬ 
icance level. This is noteworthy as different degrees of satura¬ 
tion of the C 1 41560 line are expected to produce scatter in this 
relation. This implies that the neutral-carbon content of the C i 
systems is intimately related to the reddening induced along the 
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line-of-sight or, equivalently, that the amounts of shielded gas 
and dust are tightly inter-connected. We also note that two of 
the largest three ^(B-V) values in this plot correspond to sys¬ 
tems located at Zabs ~ Zem- This may however be the result of 
small number statistics as the reddening induced by the other 
proximate systems in our sample varies substantially from one 
system to another. 

The relation between ^(B-V) and the UV bump strength, 
^bump, which we previously determined independently from 
^(B-V) (see Sect. 14.3b . is shown in the right panel of Fig. [12] 
It can be seen that both quantities are tightly correlated (b.Ocr). 
A linear least-squares fit (linear correlation coefficient of r = 
0.77), taking into account errors in both parameters, gives; ^(B- 
V) ^ 0.43 X Abump- The 2175 A extinction feature is detected 
at more than 2cr (95% confidence level) in about 30% of the 
Ci systems. In such cases, we find Abump ~ 0.4 and ^(B- 
V) ~ 0.2 mag or, equ ivalently. Ay ~ 0.6 mag. T hese values 
are comparable to what iBudzvnski & Hewetd d201 ih have found 
when targeting the strongest Mgii systems from SDSS DR 6 
[1 < Wr(/12796) < 5 A] where the 2175 A absorption is detected 
on a statistical basis only (see also Lliang et ^1201 iL for candi¬ 
date 2175 A absorption in similar systems). Interestingly, our 
measured UV bump strengths are also comparable to what has 
been observed along GRB lines-of-sight at s imilarly low leve ls 
of extinction, e.g., towards GRB 080605 (see lZafar et al.ll201^ . 

In the local Universe, 2175 A absorption can be observed 
together with reddening va lues as low as ~ 0.2 mag (see, e.g., 
iFitzpatrick & Mas^l2007h . Even at such low levels of redden¬ 
ing, the UV bump is significantly stronger along Galactic lines- 
of-sight, i.e., by up to a factor of ten, than in the present Ci 
absorber sample and/or through GRB host galaxies. This dis¬ 
crepancy may be explained if these high-redshift systems probe 
regions of the ISM affected by a star formation more vigor- 
ous than in the Gal axy. A similar argument was proposed by 
[Gordon et al.l (l2003h with the aim of explaining the variety of 
EMC and SMC extinction curves. In fact, most of the Ci ab¬ 
sorbers at the high end of the reddening tail in our sample 
are best-fit using the extinction law of the LMC2 super-shell 
near the 30 Dor star-forming regior0 (see left panel of Eig. |8]l. 
This means that the far-UV rise of the extinction curve is en¬ 
hanced and the carriers of the 2175 A absorption are depleted 
compared to Galactic lines-of-sight. This is probably the conse¬ 
quence of a high UV flux and/or the mechanical feedback from 
stars (e.g., lEischera & Dopita 1201 Ih in the vicinity of the Ci 
syste ms. In contrast, th e lack of a UV bump in typical DLAs 
(e.g.. iKhare et al.ll2012h is probably intrinsic to their low dust 
and metal contents as the lines-of-sight are likely to pass at large 
impact parameters from the absorbing galaxy. 

The main outlier in the right panel of Eig.[T2| which exhibits 
high reddening but no UV bump, is a proximate system. This 
is consistent with the above picture where the enhanced UV ra¬ 
diation field from the QSO and/or star-forming regions within 
the QSO host galaxy are expected to deplete the carriers of the 
2175 A absorption. 

7. Conclusions 

In this work, we presented a new population of QSO absorbers 
selected directly from the properties of the shielded gas, namely 

^ Super-giant shells with sizes approaching 1 kpc form the largest 
structures seen in the ISM of galaxies where large amounts of kinetic 
energy are contributed by multiple supernovae explosions and energetic 
stellar winds. 


the strongest C i absorbers, detected in low-resolution QSO spec¬ 
tra from the SDSS-II DR 7 database. These Ci absorbers, with 
lTr(T1560) > 0.4 A, are more than one hundred-times rarer than 
DLAs at Zabs = 2.5. Their number per unit redshift is increasing 
significantly below Zabs = 2, probably coupled to an increase 
in the star-formation efficiency at these redshifts. iGupta et’aH 
(1201 2h reported a similarly high detection rate of 21-cm ab¬ 
sorbers towards even lower redshifts among strong Mg ii sys¬ 
tems, which they argued must be related to the evolution of the 
CNM filling factor in the latter absorbers. 

The HI column-density distribution of C i-selected absorbers 
is flatter than that of Hi-selected absorbers. While sub-DLAs 
have much larger cross-section than DLAs, this can be under¬ 
stood as the shielding of the gas is more difficult at low Hi 
column densities and the number of clouds along the line-of- 
sight is probably smaller. Cold and dusty gas as traced by C i 
absorbers is also more likely to be found at small impact param¬ 
eters from the absorbing galaxies where a flatter V(H i) distri¬ 
bution is expected. Indeed, despite a likely bias against strong 
DLAs with large amounts of dust, we And there is among C i 
systems a probable excess of strong DLAs with log V(H i) > 21 
(atoms cm“^) compared to systematic DLA searches. This is 
reminiscent of the V(H i) distribution of DLAs within GRB host 
galaxies which is skewed t owards extremely strong DLAs (see 
flg. 10 in lLvnbo et al.ll2009ll . 

The reddening and therefore the presence of dust along the 
QSO lines-of-sight with detected C i absorption is directly re¬ 
lated to the amount of shielded gas but depends weakly on the 
total HI column density. The latter can indeed vary by more than 
a factor of ten for the same C i rest-frame equivalent width. This 
is probably the consequence of the shielded gas being clumpy 
while HI absorption samples simultaneously warm diffuse neu¬ 
tral clouds and cold, high-metallicity dusty pockets of gas. The 
presence of dust inducing significant reddening of the back¬ 
ground QSOs and/or2175 A extinction features are ubiquitous in 
about 30% of the C i absorbe rs. Several systems like these have 
been found before (see, e.g.. ISrianand et ai]|2008al: IWang et aH 
|20H . Here, we And that the UV bump is weak compared to 
Galactic lines-of-sight exhibiting the same amount of redden¬ 
ing. We interpret this as being the consequence of star formation 
in the vicinity of the systems. 

It is likely that the metal and molecular contents of Ci 
absorbers are high and actually higher than those of most 
DLAs studied till now. High-resolution spectroscopic follow¬ 
up observations of the present sample therefore opens up the 
door to systematic sea rches for carbon monoxide (CO; see 
iNoterdaeme et al.ll201 ih and molecules like CN and CH as well 
as diffuse interstellar bands at high redshift. Such a spectroscopic 
campaign will be presented in a companion paper. The typical 
reddening induced by C i absorbers along with the relation be¬ 
tween reddening and shielded-gas column density imply that the 
extinction could be high in some DLAs with C i absorption. If 
strong dusty DLAs exist, they probably have been missed in the 
current magnitude-limite d QSO samples (see also iBoisse et al.l 
ll998l:lvIadilo et a02008h . Some of the QSO lines-of-sight iden¬ 
tified here, as well as those which may be found by extending the 
present survey to even larger database^ will result in exceed¬ 
ingly long integration times on high-resolution spectrographs in¬ 
stalled on 8-lOm class telescopes. These will however be tar- 


^ Note that from the Baryon Oscillation Spectroscopic Survey, which 
is part of SDSS-III, relatively few additional C i systems are expected 
since the bulk of the new QSOs is at Zem ~ 3, which provides shorter 
C i-absorption path length. 
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gets of choice for the coming generation of Extremely Large 
Telescopes. 
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4000 4100 4200 4300 4400 4500 4600 4700 

Observed wavelength (A) 


Fig.l. SDSS spectrum of the Zem = 1-94 QSO J0815-H2640. 
The C 1 4/11560,1656 absorption lines detected at z = 1.681 are 
indicated along with several low-ionization metal lines at the 
same redshift. Among the latter, the simultaneous presence of 
Si II 41526 and Alii 41670 lines was required by our detection 
algorithm to minimize the probability of chance coincidences. In 
this system, Feii 41608 absorption is also observed even though 
relatively weak. 



WXX1656) [A] 


Fig. 2. Rest-frame equivalent widths of C i absorption lines mea¬ 
sured in SDSS spectra and associated Icr uncertainties. The up¬ 
per dashed line shows the 1:1 relation expected in case of heavily 
saturated profiles. The lower dashed line shows the expectation 
for absorption lines on the linear part of the curve-of-growth. 



W,(kl560) [A] 


Fig. 3. Completeness of our Ci search as a function of the 
C 1 41560 rest-frame equivalent width. Each point is based 
on a total of 40000 artificial Ci systems introduced in 1000 
randomly-selected SDSS QSO spectra. The vertical dotted line 
shows our chosen completeness limit of lTr(41560) = 0.4 A. 



Fig. 4. Statistical sensitivity function of our C i survey, g(z) (solid 
curve; left-hand side axis). The dips seen at z ~ 2.35 (resp. 
2.6) are due to Ci41656 (resp. Ci41560) falling over bad pix¬ 
els at the junction between the two SDSS CCD chips. Dips at 
higher redshifts (z ~ 2.8,3.05) are due to telluric lines (e.g., 
[O i] 46300). The hashed histogram (right-hand side axis) shows 
the C i-absorber number counts in different redshift bins for in¬ 
tervening systems with rest-frame equivalent width above the 
completeness limit of the survey, i.e., VTr(41560) > 0.4 A. 
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0.006 
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Fig. 5. Evolution of the number of intervening C i absorbers per 
unit redshift, iiqi, for systems with rest-frame equivalent widths 
lFr(dl560) > 0.4 A. The measurements (solid) are corrected 
for incompleteness at the latter limit. For clarity purposes, 
the uncorrected data (dashed) are displayed with a slight off¬ 
set in abscissa. The dotted curve shows the expected redshift 
behaviour of a non-evolving population. 
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A(g-r) A(r-i) A(r-z) 


Fig. 6. Upper panels: Colours of the 41 696 SDSS QSOs whose spectra were searched for Ci absorption versus QSO emission 
redshift. (g - r), (r - i) and (r - z) colours are displayed in the left, middle and right-hand side panels, respectively. The thick 
blue line shows the median QSO colour as a function of redshift. Thin blue lines are shown one standard deviation away from that 
median. The red circles indicate the 66 lines-of-sight passing through C i-bearing gas. Lower panels: Distribution of colour excesses 
[defined as, e.g.: {g - r) - {g - for A(g - r)] for the 66 QSOs whose spectra were found to exhibit C i absorption (thick red 
histogram) compared to the whole QSO sample (thin-hashed histogram, scaled down by a factor of 632 to have the same area). In 
each panel [from left to right: A(g - r), A(r - i) and A(r - z)], vertical dotted red lines indicate the median colour excesses of the 
sample of QSOs with detected C i absorption. The median values for the whole QSO sample are shown by vertical dashed lines. 
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1500 


Rest-frame wavelength (A) at 1.681 
2000 2500 3000 




Fig. 7. Left panel: Illustration of reddening meas urement based on OSO S EP fitting. The blue curve shows our best fit to QSO 
continuum using the QSO template spectrum from IVanden Berk et all (1200 Ih and the best-fitting LMC extinction law at the redshift 
of the CI absorber. The non-extinguished QSO template is displayed in grey. The additional presence of a 2175 A absorption feature 
is shown by the light-red shaded area (i.e., Abump)- The broken orange line at zero ordinate indicates the wavelength regions devoid 
of QSO emission lines, strong absorption lines and sky residuals, which were used during the minimization. Upper right panel: 
Distribution of reddening measured in the QSO control sample assuming either a SMC (red) or the best-fitting extinction law (in 
this case, LMC; blue). The reddening values of the sight-line under consideration are indicated by downward arrows. Lower right 
panel: Same as above but for the distribution of 2175 A bump strengths. 




E(B-V) [mag] 


Fig. 8 . Left panel: Histogram of reddening for the sample of 66 QSOs whose spectra were found to exhibit C i absorption (multi¬ 
color). The colour/hash coding for the different C i-detected lines-of-sight relates to the best-fitting extinction law (red; SMC; blue: 
LMC; green; LMC2; purple: MW). The distribution of reddening from the QSO control sample (C.S.), i.e., the sum of the normalized 
distributions of individual QSO control samples, is displayed with dashed lines. Right panel: Same as in the left panel but for the 
histogram of 2175 A bump strengths. The distribution from the QSO control sample is well represented by a Gaussian profile (solid 
curve). 
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log N(HI) [atoms cm 

Fig. 9. Neutral atomic-hydrogen column-density distribution 
of Ci absorbers (red-filled histogram; left-hand side axis). 
Proximate systems are displayed in salmon. The N(H i) distri¬ 
bution of intervening DLAs from SDSS DR7 (iNoterdaeme et al.l 
l2009bl blue-hashed histogram) is over-plotted using a different 
scaling (right-hand side axis) so that the areas of the two his¬ 
tograms above log N(H i) = 20.3 are equal. The blue dashed- 
dotte d curve is a fit o f the d istribution in the sub-DLA regime 
from iProchaska et al.l (l2014h . The inset shows the cumulative 
A(Hi) distributions of the intervening Ci and DLA samples 
starting at logA(Hi) = 20.3 (atoms cm“^) and the two-sided 
Kolmogorov-Smirnov test probability that the two distributions 
come from the same parent population. 



Fig. 10. C 1/11560 rest-frame equivalent width versus A(Hi), 
the neutral atomic-hydrogen column density of Ci absorbers. 
Intervening (resp. proximate) systems are displayed as filled 
(resp. empty) black circles. Data points are encircled in red as 
they correspond to the sub-sample of C i systems with measured 
A(H i), i.e., the H i sub-sample. The dashed (resp. dashed-dotted) 
curve materializes the limit above which systems are expected 
to have in excess of one-tenth of Solar (resp. in excess of Solar) 
metallicity (see text). The completeness limit of the survey is 
indicated by the horizontal line. 
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WXX1560) [A] 



-0.4 -0.2 


^bump 


Fig. 12. Left panel: Reddening of the background QSOs with detected C i absorption (overall sample) versus C i /il560 rest-frame 
equivalent width. Symbol conventions are the same as in Fig. [10] Data points encircled in red correspond to the H i sub-sample. 
Right panel: Same as in the left panel with the exception that reddening is plotted against 2175 A bump strength. 



Fig. 11. Reddening of the background QSOs with detected C i 
absorption versus A(Hi) of the Ci absorbers. Symbol conven¬ 
tions are the same as in Fig.[T0| The solid line shows the relation 
observed in the loc al ISM where £(B-V)/A (H) = 1.63 x 10 ^2 
mag atoms * cm2 dGudennavar et al.ll2012h . The dotted line il¬ 
lustrates a reddening per hydrogen atom ten times higher than 
that. The dashed cu rves correspond to the observations of the 
Magellanic Clouds (iGordon et alJl2003l) . The long-dashed line 
shows the relation derived for typical DLAs at z^bs ~ 2.8 (e.g., 
IVladilo et al.ll2()()i l^ are et IDI 20 T 1 . 
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Table 1. C i systems identified in SDSS-DR7 QSO spectra. 


QSO name 

plate 

MJD 

fibre 

■^em 

'^abs 

lV,('il560)“ 

[A] 

lV,(dl656)“ 

[A] 

S/N*’ 

logA(Hl)^,Dss 

[atoms cm 

logA(Hl)-vEs 

[atoms cm 

J021606.13-002104.9 

0405 

51816 

114 

2.22 

1.737 

0.33 ±0.10 

0.77 ± 0.08 

14 



J 030049.24-072137.8 

0458 

51929 

437 

2.11 

1.536 

0.79 ±0.15 

1.17 ±0.21 

5 



J080801.74-H330009.2 

0825 

52289 

166 

1.90 

1.888 

0.24 ± 0.05 

0.39 ± 0.05 

18 



J081116.06-H083837.7 

2570 

54081 

566 

2.07 

1.906 

0.23 ± 0.06 

0.21 ± 0.06 

15 


20.10 ±0.10 

J081540.60-t-264021.6 

1266 

52709 

355 

1.94 

1.681 

1.70 ±0.21 

1.33 ±0.16 

5 



J 082003.404-155932.9 

2272 

53713 

378 

1.95 

1.547 

0.72 ±0.15 

0.77 ±0.11 

7 



J085206.65-H 193548.4 

2281 

53711 

452 

1.91 

1.788 

0.34 ± 0.06 

0.86 ± 0.07 

13 


20.00 ± 0.20 

J085437.60-H031734.8 

0565 

52225 

480 

2.24 

1.567 

0.40 ± 0.09 

0.43 ± 0.07 

12 



J085726.79-t-185524.3 

2281 

53711 

630 

1.89 

1.730 

0.20 ± 0.02 

0.43 ± 0.04 

33 


19.60 ± 0.20 

J090558.75-t-553039.2 

0450 

51908 

470 

1.57 

1.579 

0.13 ±0.05 

0.15 ±0.05 

16 



J090942.56-H532126.5 

0553 

51999 

425 

2.07 

1.628 

0.47 ± 0.08 

0.70 ± 0.09 

12 



J091516.27-H071750.9 

1194 

52703 

506 

2.38 

1.681 

0.50 ±0.10 

0.95 ±0.13 

11 



J091721.37-f015448.1 

0473 

51929 

426 

2.18 

2.107 

0.67 ±0.16 

0.93 ±0.16 

5 



1092759.794-154321.8 

2579 

54068 

364 

1.81 

1.731 

0.61 ±0.10 

0.69 ± 0.09 

9 


21.10 ±0.20 

J 094529.37-4552525.7 

0556 

51991 

045 

2.24 

1.867 

0.74 ±0.12 

0.91 ±0.13 

7 



J 095114.97-4545736.5 

0769 

52282 

458 

1.80 

1.613 

0.77 ± 0.07 

0.78 ± 0.06 

14 



J 101544.34-4310617.2 

1953 

53358 

604 

1.56 

1.596 

0.20 ± 0.09 

0.34 ± 0.09 

9 



J 104705.76-4205734.5 

2478 

54097 

558 

2.01 

1.775 

0.59 ±0.14 

0.96 ±0.13 

7 



J 105436.96-4542343.7 

0907 

52373 

129 

1.50 

1.513 

0.49 ± 0.08 

0.87 ± 0.08 

10 



J 105746.42-4662250.5 

0490 

51929 

179 

2.34 

2.088 

0.30 ± 0.06 

0.48 ± 0.06 

13 



J 111756.53-4143716.9 

1753 

53383 

353 

2.09 

2.001 

0.67 ± 0.09 

0.65 ± 0.08 

10 


19.80 ±0.15 

J 112250.95-4143732.5 

1753 

53383 

477 

2.05 

1.554 

0.17 ±0.04 

0.19 ±0.03 

28 



J 112913.29-023740.9 

0327 

52294 

308 

1.86 

1.623 

0.60 ±0.10 

0.99 ±0.12 

8 



J 113341.29-005740.1 

0282 

51658 

215 

1.68 

1.706 

1.36 ±0.25 

1.24 ±0.24 

4 



J 114141.84-4444206.1 

1367 

53083 

499 

1.96 

1.903 

0.48 ±0.11 

0.76 ±0.13 

7 



J 115153.87-4150945.0 

1762 

53415 

495 

3.05 

2.400 

0.36 ±0.11 

0.37 ±0.11 

7 

20.15 ±0.33= 

20.15 ±0.10 

J 115705.51-4615521.7 

0777 

52320 

107 

2.51 

2.460 

1.22 ±0.14 

1.30 ±0.14 

7 

21.80 ±0.20* 


J 120935.79-4671715.7 

0493 

51957 

178 

2.03 

1.843 

0.47 ±0.11 

0.48 ±0.11 

8 



J 122825.67-4303038.6 

2234 

53823 

430 

3.91 

3.098 

0.17 ±0.07 

0.26 ± 0.07 

11 

21.05 ±0.21= 


J 123714.61-4064759.6 

1628 

53474 

193 

2.78 

2.691 

0.37 ± 0.06 

0.77 ± 0.08 

14 

20.15 ±0.28= 

20.00 ±0.158 

J 124708.43-4500320.8 

1278 

52735 

039 

2.27 

2.135 

1.05 ±0.10 

1.42 ±0.11 

9 



J 124829.71-4284858.1 

2239 

53726 

219 

1.54 

1.513 

0.32 ± 0.07 

0.65 ± 0.07 

12 



J 124841.45-4302433.0 

2239 

53726 

469 

2.06 

1.691 

0.57 ± 0.06 

0.87 ± 0.06 

15 



J 125552.60-4223424.4 

2649 

54212 

507 

1.82 

1.526 

0.25 ±0.11 

1.18 ±0.28 

8 



J 130225.28-4211158.6 

2650 

54505 

296 

1.76 

1.656 

0.56 ±0.12 

0.49 ±0.10 

8 



J 130628.87-4281550.8 

2009 

53904 

288 

2.10 

2.012 

0.28 ±0.10 

0.56 ±0.10 

8 



J 130828.43-4584000.6 

0958 

52410 

224 

3.09 

2.473 

0.39 ± 0.07 

0.61 ±0.14 

11 

20.00 ±0.15 


J 131129.11-4222552.6 

2651 

54507 

398 

3.14 

3.093 

0.33 ± 0.05 

0.51 ±0.08 

16 

20.50 ±0.10 

20.62 ±0.10 

J 131400.57-4054319.5 

0850 

52338 

620 

1.89 

1.583 

0.23 ±0.10 

0.70 ±0.10 

8 



J 134122.51-4185214.0 

2642 

54232 

494 

2.00 

1.544 

0.16 ±0.03 

0.15 ±0.04 

27 



J 134601.10-4064408.4 

1803 

54152 

074 

2.09 

1.512 

0.79 ± 0.20 

1.33 ±0.28 

6 



J 135122.00-4461828.5 

1466 

53083 

416 

1.81 

1.606 

0.99 ±0.15 

1.58 ±0.17 

10 



J 141550.47-4300146.9 

2129 

54252 

288 

2.08 

1.676 

1.35 ±0.10 

1.61 ±0.12 

13 



J 141606.79-4180403.2 

2759 

54534 

359 

2.13 

1.622 

0.47 ±0.11 

0.80 ±0.13 

8 



J 143243.93-4330746.7 

1841 

53491 

513 

2.09 

2.058 

0.28 ± 0.05 

0.44 ± 0.05 

17 



J 143657.87-4291100.6 

2138 

53757 

259 

1.77 

1.596 

0.35 ±0.11 

0.56 ±0.11 

8 



J 143912.05-4111740.6 

1711 

53535 

374 

2.58 

2.418 

0.50 ± 0.08 

0.68 ± 0.09 

10 

20.27 ± 0.24= 

20.10 ±0.10*’ 

J 144929.27-4333811.0 

1646 

53498 

065 

2.17 

2.021 

0.47 ±0.12 

0.62 ±0.12 

7 



J 145432.54-4343523.9 

1384 

53121 

269 

1.61 

1.580 

0.59 ±0.16 

0.59 ±0.13 

6 



J 145953.25-4012944.2 

0538 

52029 

039 

1.66 

1.623 

0.36 ±0.17 

0.62 ±0.16 

5 



J 150738.73-4415530.6 

1291 

52735 

505 

1.79 

1.674 

0.50 ±0.11 

0.85 ±0.16 

7 



J 152209.12-4083020.0 

1721 

53857 

192 

1.93 

1.627 

0.32 ±0.10 

0.71 ±0.11 

8 



J 160320.76-4170117.8 

2200 

53875 

318 

1.99 

1.890 

0.16 ±0.07 

0.43 ± 0.08 

12 


19.20 ±0.15 

J 160457.51-4220300.5 

2205 

53793 

328 

1.98 

1.641 

1.34 ±0.12 

1.60 ±0.13 

7 



J 161526.65-4264813.8 

1576 

53496 

394 

2.18 

2.118 

0.28 ± 0.06 

0.56 ±0.12 

15 


20.55 ±0.15 

J 162321.43-4135532.4 

2202 

53566 

132 

1.75 

1.751 

0.39 ± 0.09 

0.60 ±0.10 

9 


19.95 ±0.15 

J 164610.20-4232923.0 

1424 

52912 

301 

2.06 

1.998 

0.36 ± 0.05 

0.58 ± 0.08 

15 


19.68 ±0.15 

J 170542.92-4354340.4 

0974 

52427 

561 

2.01 

2.038 

0.70 ± 0.08 

1.40 ±0.10 

10 



J 212329.47-005053.0 

0987 

52523 

103 

2.26 

2.060 

0.11 ±0.02 

0.11 ±0.02 

41 


19.25 ± 0.20‘ 


“ Rest-frame equivalent width of the absorption detected at the position of the C i lines (measurements drawn from SDSS spectra). 
Average signal-to-noise ratio per pixel around the C i lines. 

Neutral atomic-hydrogen column density of the systems measured from Lya absorption in SDSS spectra. 

A(H l) values derived from follow-up UVES sp ec troscopy. _ 

'^iNoterdaeme et al] d2009bl) : ^ IWang et'H] (l2012h ; ^ iNoteriaeme et alj(l2010al) : ‘‘ iNoterdaeme et al.1 d2008d) : ‘ iMilutinovic et alj(l2010h : t 
iLedoux et al.l 1200^. 
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Table 1. continued. 


QSO name 

plate 

MJD 

fibre 

■Sem 

2abs 

W'r(dl560)“ 

[A] 

W',(^1656)“ 

[A] 

S/N*’ 

logA(Hl)^,j,,3 
[atoms cm 

logA(Hl)-^E, 
[atoms cm 

J222910.16-t-141402.2 

0738 

52521 

382 

2.11 

1.586 

0.33 ±0.13 

0.47 ±0.13 

6 



J 225719.04-100104.7 

0724 

52254 

150 

2.08 

1.836 

0.64 ± 0.05 

0.89 ± 0.05 

17 


20.38 ± 0.07 

J 233133.05-090246.6 

0646 

52523 

619 

2.44 

1.734 

0.41 ± 0.04 

0.35 ± 0.04 

20 



J 233156.49-090802.0 

0646 

52523 

616 

2.66 

2.143 

1.34 ±0.13 

1.17±0.11 

8 


21.15 ±0.15 

J233633.81-105841.5 

0647 

52553 

201 

2.04 

1.829 

0.22 ± 0.05 

0.40 ± 0.07 

15 



J 234023.67-005327.1 

0385 

51877 

204 

2.09 

2.054 

0.11 ±0.02 

0.17 ±0.02 

37 


20.33 ± 0.03 

J 235057.87-005210.0 

0386 

51788 

137 

3.03 

2.426 

0.62 ± 0.08 

0.63 ± 0.08 

11 

20.39 ± 0.27= 

20.50 ±0.10i 


“ Rest-frame equivalent width of the absorption detected at the position of the C i lines (measurements drawn from SDSS spectra). 
Average signal-to-noise ratio per pixel around the C i lines. 

Neutral atomic-hydrogen column density of the systems measured from Lya absorption in SDSS spectra. 

A(H l) value derived from f ollow-up UVES spe ct roscopy. ___ 

'^iNoterdaeme et al] d2009bl) : ^ IWang (1201^ ; ^ iNoterdaeme et alj(l2010ah : ‘‘ iNoterdaeme et al.1 d2008d) : ‘ iMilutinovic et alj(l2010h : j 
iLedoux et al.l 1200^. 

Table 2. Intervening C i-absorber number counts. 


CI samples 

Wr(dl560) 

[A] 

2min 

^max 

Az 

#of Ci 

4 1 

(xlO-3) 

All systems 

> 0.4 

1.5 

1.9 

6268 

22 

3.5 ±0.8 



1.9 

3.35 

6710 

8 

1.2 ±0.4 

Strongest 

> 0.64 

1.5 

1.9 

6268 

10 

1.6 ±0.5 



1.9 

3.35 

6710 

4 

0.6 ± 0.3 

Weaker 

> 0.4, < 0.64 

1.5 

1.9 

6268 

12 

1.9 ±0.6 



1.9 

3.35 

6710 

4 

0.6 ± 0.3 


“ Prior to correcting for incompleteness (i.e., 18% at M4(/11560) = 0.4 A). 
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Table 3. SED fitting of QSO spectra with detected C i absorbers. 


QSO name 

■^em 

^abs 

£(B-V) 

[mag] 

best fit 

^bump 

C.S.^: 

Nqso 

E(B-Y) [mag] 
median cr 

^bump 

median cr 

J 021606.13-002104.9 

2.22 

1.737 

0.039 

SMC 

0.097 

173 

-0.008 

0.053 

-0.009 

0.123 

J 030049.24-072137.8 

2.11 

1.536 

0.075 

SMC 

0.160 

448 

-0.004 

0.048 

0.002 

0.136 

J 080801.74-t330009.2 

1.90 

1.888 

0.029 

SMC 

-0.032 

151 

0.001 

0.044 

-0.028 

0.112 

J081116.06-H083837.7 

2.07 

1.906 

0.081 

SMC 

-0.027 

247 

-0.003 

0.049 

-0.031 

0.127 

J081540.60-t-264021.6 

1.94 

1.681 

0.109 

LMC 

0.364 

428 

-0.012 

0.043 

-0.027 

0.142 

J082003.40-H 155932.9 

1.95 

1.547 

0.203 

LMC2 

0.420 

373 

-0.004 

0.069 

-0.012 

0.142 

J085206.65-H 193548.4 

1.91 

1.788 

0.050 

SMC 

0.032 

580 

-0.003 

0.044 

-0.020 

0.112 

J085437.60-H031734.8 

2.24 

1.567 

0.070 

SMC 

0.263 

192 

-0.003 

0.054 

0.050 

0.106 

J085726.79-tl85524.3 

1.89 

1.730 

0.016 

SMC 

0.063 

81 

-0.004 

0.053 

-0.064 

0.107 

J 090558.754-553039.2 

1.57 

1.579 

-0.019 

SMC 

0.070 

692 

-0.002 

0.041 

0.011 

0.080 

J090942.56-t532126.5 

2.07 

1.628 

0.003 

SMC 

0.003 

527 

-0.005 

0.044 

-0.019 

0.140 

J091516.27-H071750.9 

2.38 

1.681 

-0.010 

SMC 

0.117 

40 

-0.017 

0.077 

0.019 

0.206 

J091721.37-t-015448.1 

2.18 

2.107 

0.121 

SMC 

0.227 

86 

-0.014 

0.031 

0.006 

0.114 

J092759.79-tl54321.8 

1.81 

1.731 

0.217 

LMC2 

0.331 

357 

-0.008 

0.059 

-0.018 

0.092 

J094529.37-H552525.7 

2.24 

1.867 

0.120 

SMC 

0.015 

314 

-0.006 

0.043 

-0.016 

0.130 

J095114.97-H545736.5 

1.80 

1.613 

-0.004 

SMC 

-0.027 

397 

-0.004 

0.042 

-0.015 

0.109 

J 101544.344-310617.2 

1.56 

1.596 

0.016 

SMC 

-0.006 

685 

0.000 

0.043 

0.012 

0.085 

J 104705.764-205734.5 

2.01 

1.775 

0.159 

LMC2 

0.277 

282 

-0.021 

0.055 

-0.036 

0.133 

J 105436.96-H542343.7 

1.50 

1.513 

0.014 

SMC 

-0.008 

319 

0.001 

0.037 

0.007 

0.072 

J 105746.42-H662250.5 

2.34 

2.088 

0.040 

SMC 

0.112 

277 

-0.010 

0.037 

-0.046 

0.120 

J 111756.534-143716.9 

2.09 

2.001 

0.049 

SMC 

0.085 

560 

-0.005 

0.036 

-0.030 

0.097 

J 112250.954-143732.5 

2.05 

1.554 

-0.028 

SMC 

0.032 

106 

0.003 

0.047 

-0.026 

0.138 

J 112913.29-023740.9 

1.86 

1.623 

0.020 

SMC 

-0.151 

603 

-0.003 

0.045 

-0.026 

0.124 

J 113341.29-005740.1 

1.68 

1.706 

0.266 

SMC 

0.062 

386 

-0.002 

0.042 

0.002 

0.094 

J114141.84-t444206.1 

1.96 

1.903 

0.149 

LMC2 

0.368 

398 

-0.023 

0.060 

-0.028 

0.109 

J 115153.874-150945.0 

3.05 

2.400 

-0.007 

SMC 

0.110 

92 

-0.007 

0.033 

-0.105 

0.176 

J 115705.514-615521.7 

2.51 

2.460 

0.215 

SMC 

0.370 

148 

-0.014 

0.033 

-0.087 

0.101 

J 120935.79-4671715.7 

2.03 

1.843 

0.094 

MW 

0.394 

300 

-0.013 

0.028 

-0.039 

0.119 

J 122825.67-4303038.6 

3.91 

3.098 

-0.034 

SMC 

0.067 

43 

-0.008 

0.061 

0.412 

0.471 

J 123714.614-064759.6 

2.78 

2.691 

0.143 

LMC2 

0.431 

35 

-0.026 

0.074 

-0.168 

0.156 

J 124708.434-500320.8 

2.27 

2.135 

0.075 

LMC 

0.319 

357 

-0.012 

0.029 

-0.037 

0.124 

J 124829.71-4284858.1 

1.54 

1.513 

-0.035 

SMC 

-0.049 

1035 

-0.002 

0.041 

0.022 

0.092 

J 124841.45-4302433.0 

2.06 

1.691 

-0.025 

SMC 

-0.019 

583 

-0.004 

0.042 

-0.024 

0.136 

J 125552.604-223424.4 

1.82 

1.526 

0.036 

SMC 

0.141 

898 

-0.004 

0.041 

-0.015 

0.134 

J 130225.284-211158.6 

1.76 

1.656 

-0.004 

SMC 

-0.139 

452 

-0.012 

0.036 

-0.002 

0.114 

J 130628.87-4281550.8 

2.10 

2.012 

0.078 

SMC 

0.015 

465 

-0.008 

0.034 

-0.030 

0.099 

J 130828.43-4584000.6 

3.09 

2.473 

0.020 

SMC 

-0.158 

54 

-0.010 

0.037 

-0.144 

0.211 

J 131129.11-4222552.6 

3.14 

3.093 

0.036 

SMC 

0.067 

82 

-0.011 

0.036 

-0.054 

0.345 

J 131400.574-054319.5 

1.89 

1.583 

0.046 

SMC 

0.270 

340 

-0.006 

0.044 

-0.013 

0.120 

J 134122.51-4185214.0 

2.00 

1.544 

-0.016 

SMC 

0.000 

32 

0.003 

0.053 

-0.001 

0.185 

J 134601.10-4064408.4 

2.09 

1.512 

0.015 

SMC 

0.034 

118 

-0.015 

0.040 

0.036 

0.130 

J 135122.00-4461828.5 

1.81 

1.606 

0.142 

LMC2 

0.408 

680 

-0.005 

0.067 

-0.027 

0.123 

J 141550.47-4300146.9 

2.08 

1.676 

0.074 

SMC 

0.226 

83 

-0.001 

0.045 

-0.033 

0.245 

J 141606.79-4180403.2 

2.13 

1.622 

0.056 

SMC 

-0.027 

108 

-0.016 

0.035 

0.018 

0.150 

J 143243.93-4330746.7 

2.09 

2.058 

0.010 

SMC 

-0.005 

352 

-0.003 

0.041 

-0.039 

0.095 

J 143657.874-291100.6 

1.77 

1.596 

0.147 

LMC 

0.514 

812 

-0.003 

0.049 

-0.014 

0.111 

J 143912.054-111740.6 

2.58 

2.418 

0.172 

LMC2 

0.271 

52 

-0.033 

0.145 

-0.080 

0.098 

J 144929.27-4333811.0 

2.17 

2.021 

0.018 

SMC 

-0.157 

89 

-0.009 

0.029 

0.025 

0.137 

J 145432.54-4343523.9 

1.61 

1.580 

0.293 

LMC2 

0.631 

346 

-0.002 

0.058 

-0.002 

0.089 

J 145953.25-4012944.2 

1.66 

1.623 

0.157 

LMC2 

0.331 

842 

-0.006 

0.057 

0.011 

0.089 

J 150738.73-4415530.6 

1.79 

1.674 

0.056 

LMC 

0.251 

806 

-0.005 

0.040 

-0.003 

0.111 

J 152209.12-4083020.0 

1.93 

1.627 

0.017 

SMC 

-0.057 

169 

-0.016 

0.035 

-0.031 

0.134 

J 160320.76-4170117.8 

1.99 

1.890 

-0.010 

SMC 

0.030 

406 

-0.004 

0.043 

-0.033 

0.108 

J 160457.514-220300.5 

1.98 

1.641 

0.263 

LMC2 

0.735 

225 

-0.007 

0.075 

-0.041 

0.151 

J 161526.654-264813.8 

2.18 

2.118 

0.035 

SMC 

0.120 

324 

-0.005 

0.042 

-0.037 

0.095 

J 162321.43-4135532.4 

1.75 

1.751 

0.031 

SMC 

0.034 

816 

-0.005 

0.036 

0.006 

0.093 

J 164610.20-4232923.0 

2.06 

1.998 

-0.027 

SMC 

-0.031 

284 

-0.004 

0.043 

-0.037 

0.092 

J 170542.924-354340.4 

2.01 

2.038 

0.146 

LMC2 

0.364 

210 

-0.009 

0.075 

-0.030 

0.077 

J 212329.47-005053.0 

2.26 

2.060 

-0.011 

SMC 

0.000 

32 

0.002 

0.041 

-0.023 

0.123 

J222910.16-4141402.2 

2.11 

1.586 

-0.024 

SMC 

0.026 

538 

-0.004 

0.041 

0.006 

0.136 

J225719.04-100104.7 

2.08 

1.836 

0.045 

SMC 

0.177 

192 

-0.007 

0.039 

-0.040 

0.129 

J 233133.05-090246.6 

2.44 

1.734 

0.101 

LMC 

0.393 

38 

0.002 

0.089 

0.084 

0.215 

J 233156.49-090802.0 

2.66 

2.143 

0.136 

SMC 

0.234 

46 

-0.008 

0.078 

-0.015 

0.104 

J233633.81-105841.5 

2.04 

1.829 

0.004 

SMC 

-0.103 

484 

-0.003 

0.045 

-0.026 

0.123 

J 234023.67-005327.1 

2.09 

2.054 

0.057 

SMC 

-0.023 

34 

-0.001 

0.047 

-0.071 

0.113 

J 235057.87-005210.0 

3.03 

2.426 

0.003 

SMC 

0.001 

42 

-0.013 

0.067 

-0.098 

0.169 


C.S. = control sample. 







